Longitudinal Laser Shaping in Laser Wakefield Accelerators 
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response can be written in terms of the independent di- 
mensionless variables C = kp{z — Vgt) and r = kpct^ 
where kp is the plasma wavenumber, and Vg ~ —cis the 
laser group velocity (for convenience the laser is mov- 
ing in the negative z direction). Introducing dimension- 
less normalized scalar and vector potentials ipiO ^'^d 
a(C), the parallel and perpendicular electric fields are 
= —{m(?kple)d(i>ldQ and = —{mc/e)da/dt — 
— {mc^kp/e)da/dC,. The wakefield generation equation 
can then be written as 



i/i + fl'(C) 
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We study the energetics of wake excitation during the laser-plasma interaction in application to 
laser wakefield accelerators. We find that both the wake amplitude and the accelerating efficiency 
(transformer ratio) can be maximized by properly shaping the longitudinal profile of the driving laser 
pulse. The corresponding family of laser pulse shapes is derived in the nonlinear regime of laser- 
plasma interaction. Such shapes provide theoretical upper limit on the magnitude of the wakefield 
and efficiency by allowing for uniform photon deceleration inside the laser pulse. We also construct 
realistic optimal pulse shapes that can be produced in finite-bandwidth laser systems. 
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Current plasma accelerators can be broadly subdivided 
into two classes - the laser driven scheme [Q, and the 
electron beam driven scheme While the accelerat- 
ing plasma waves excited in these two schemes are very 
similar, the physics of wakefield excitation is quite differ- 
ent between them. In the beam-driven Plasma Wakefield 
Accelerator (PWFA) the electron beam loses energy to 
the plasma through interaction with the induced elec- 
trostatic field, while in the laser-driven Laser Wakefield 
Accelerator (LWFA) energy loss occurs via photon red- 
shift or deceleration This process is due to the pres- 
ence of wake-induced stationary modulations in the re- 
fractive index of the plasma as seen in the laser comoving 
frame Q. With the recent advances in laser technology, 
the creation of ultrashort pulses with virtually arbitrary 
temporal profiles is now possible by using amplitude and 
phase pulse shapers . Since perturbations to the refrac- 
tive index of the plasma depend on pulse shape, different 
laser shapes will vary in their coupling to the plasma. A 
natural question therefore arises: what laser shape is the 
"optimal" for laser-plasma accelerators? 

The number of parameters involved in selecting a par- 
ticular pulse shape can be overwhelming. One can char- 
acterize a shape by the value of its total energy, length, 
maximum laser field, amplitude of plasma wake, etc., in 
addition to an infinite-dimensional space of actual shape 
functions. Luckily, not all of these parameters are inde- 
pendent or even well-defined. In this Letter we argue that 
the only two meaningful parameters that describe a laser 
shape from the stand point of wake excitation are the 
total pulse energy and its depletion length. Using these 
parameters different laser shapes can be consistently clas- 
sified and cross-compared while desired properties such 
as wake amplitude or efficiency can be optimized. 

Let us consider a homogeneous unmagnetized plasma 
which is charge neutral, at rest, and has initial density 
Up in the absence of electromagnetic wave. Laser propa- 
gates along the z axis with initial frequency loq ujp = 
yjAmpTLpJnif,. In the laser comoving frame, the plasma 
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where rif, is the perturbed electron plasma density, x = 
1 -f is the modified electrostatic potential, and a^(C) 
is the dimensionless laser intensity envelope averaged 
over fast oscillations. Prior to the arrival of the laser 
the normalized wakefield £ = eE\\/meCLOp — —dx/dC, 
is zero. A formal solution for the electric field out- 
side the laser can be written as the first integral of (|^): 
[£°^\C)Y = -{x - l)^/x+ a^x'/x^dC, which reaches 
a maximum value at x = 1: 
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This expression can be understood in terms of the de- 
position of laser energy into plasma. Due to negligible 
scattering, the photon number in the laser is essentially 
constant, and the laser energy loss is predominantly in 
the form of frequency redshift, or photon deceleration 



duj I ujp ^ d rie k ( ^ 

dz 2 Lu d( ^rip 2lo ^\dCx 



(3) 



The energy density in the wake from (j^) can then be 
interpreted as the intensity-weighted integral of the pho- 
ton deceleration throughout the pulse. Let's denote the 
wake-dependent part of the photon deceleration function 
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as k{() = x' jx^ . It is closely related to the character- 
istic laser depletion length ld^ or the distance in which 
the maximally decelerated laser slice redshifts down to 

assuming a nonevolving wakefield. From (^), 1^ = 
[{luq/uJp)'^ — l]/kpKmax, where Kmax is the maximum of 
inside the laser. The value of the peak wakefield in 

is bounded from above by the total laser energy and 
the maximum photon deceleration: 
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One possible optimization problem can then be formu- 
lated as follows: given some fixed laser energy Eq (i-S-, 
the integral of a^), and the desired depletion length (i.e., 
Kmax = kq), what lascr shape would produce the largest 
possible wakefield? From (Q) it is clear that such a shape 
should maintain a constant maximal photon deceleration 
throughout the pulse. If the laser is present for C > Oj 
then in order to satisfy the boundary conditions the pho- 
ton deceleration should rise from to value kq at the 
very beginning of the pulse, e.g., like a step- function: 
k(C) = Ko6{('^). Here, = C — 0+ in order to avoid 
ambiguities with the values of step- function at 0. The 
corresponding laser profile is then found from the wake 
equation (0): 
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where C £ [0, C/ < 1/ko], and S{(+) 
rC>o 



is a delta-function 



such that J^^" S{y^)dy = 1. A schematic drawing of 
the optimal laser intensity variation and its associated 
plasma wakefield are shown in Fig. |l|. Generally, the 
shape consists of a (5-function at the front followed by a 
ramp in intensity which is cut off at Q. In the linear 
regime, when ^ 1, kq — > 0, the ramp reduces to a tri- 
angular shape found in §,0: = 2ko((5(C+) + C)- We 
note that (||) describes a family of shapes, rather than 
a fixed shape. The actual profile of the optimal pulse 
depends on the deceleration parameter kq set by the de- 
sired depletion length and the pulse length Q, which is 
determined from the available total energy: 
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Although the pulse length cannot exceed = l/'*Oi the 
rise of towards the end of the pulse guarantees that 
any finite laser energy can be accommodated for Q < Cc. 
The two terms in represent the energy contained in 
the (5-function precursor and the main pulse. It is clear 
that for a fixed total energy there exists a maximum value 
of Kq = eo/2 which is achieved when Q —> 0, i.e., all 
of the energy is concentrated in the (5-function. This 
shape, which is a particular case of the general optimal 
shape (||), excites the largest possible wakefield and has 



the smallest depletion length among all pulses of fixed 
energy. For circularly polarized pulses with cylindrical 
transverse crossection of radius ro and wavelength A, the 
maximum achievable wake is then given by: 
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where C/q is the total pulse energy (in Joules) and E^i, = 
96[np/10 
ing field. 



96[np/10^*cm •^]GV/m is the nonrelativistic wavebreak 




FIG. 1. General shape of the nonlinear optimal laser inten- 
sity profile and its corresponding wakefield (arbitrary units). 

The method used for obtaining the optimal shape (|^) 
is actually more general and can be used to determine 
laser shapes that generate other variations in the nonlin- 
ear index of refraction. Having a physical requirement 
for the refractive index r] = [1 — {ujp/uj)'^ne/^np]^/'^ = 
[1 — {uip/ijoYl/xY/"^, which for this case is the require- 
ment of uniformity of photon deceleration, provides a 
constraint on a;(C), which can then be used to find the 
laser shape from the wake equation. Alas, such a "re- 
verse" solution is not always guaranteed to yield a physi- 
cal (i.e., positive) a^(C), so, in general, caution is advised. 

While the generation of large accelerating gradients is 
a prerequisite for a successful accelerating scheme, the ef- 
ficiency of acceleration should also be considered. For an 
accelerating scheme that involves transfer of energy from 
the driver beam to the accelerating beam, efficiency is 
measured in terms of the transformer ratio, or the ra- 
tio of the maximum rate of energy gain per particle of 
accelerating beam to the maximum rate of energy loss 
per particle of the driving beam. In the case of laser- 
plasma accelerators, where the driving and accelerating 
beams consist of particles of different species, the follow- 
ing kinematic definition is more useful: 
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where 7^ and 7^ are Lorentz factors for the driving and 
accelerating beams. In LWFA the particles in the trailing 
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electron bunch are accelerated via electrostatic interac- 
tion with the wake, so \d^a/dz\max = \eE^°-'^\/mf,c^ = 
^pl^maxl- F'o'^ laser propagating in plasma 7^ « uj/ujp, 
so \d^d/dz\ is the photon frequency shift given by (0). 
The transformer ratio for LWFA is then: 
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It follows from this definition that the transformer ratio 
is not only a measure of local accelerating efficiency, but 
is also directly related to the maximum energy that can 
be transferred to the accelerating beam particle over the 
laser depletion length (assuming no evolution of the wake 
during propagation). 

As there is no absolute maximum of R, we can only 
look for a shape that maximizes R subject to constraints. 
For instance, among the pulses of fixed energy and deple- 
tion length R is maximized by a pulse that produces the 
largest wakefield. But this is precisely the optimal shape 
found above. A more general proof involves considering 
all pulses irrespective of total energy that create a given 
value of wakefield. It can then be shown that a pulse 
which has the largest depletion length among these must 
maintain a constant photon deceleration inside the pulse, 
which again points to the shape (||). From (||) the opti- 
mal transformer ratio is then: 
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where Lp — C,f/kp is the pulse length. In the linear 
regime optimal transformer ratios for both LWFA and 
PWFA schemes scale identically with the pulse/beam 
le ngth: fl^WFA ^ {2lo /ujp)^l + {kpLpuUe?, i?^^^^ ^ 
\/\ + {kpLbeamY • The LWFA scheme is intrinsically 
more efficient by a factor of 2uj /ujp^ which is needed for vi- 
ability of LWFA since lasers are typically "slower" drivers 
than electron beams. 

The advantage of using the optimal pulse shape is best 
seen in comparison with the unshaped (Gaussian) pulse. 
For a given Gaussian pulse (or any other non-optimal 
shape) one can always construct a corresponding opti- 
mally shaped pulse with the same laser energy such that 
the photon deceleration across the optimal pulse equals 
to the peak photon deceleration in the unshaped one (i.e., 
both pulses have equal depletion lengths). Unshaped 
pulses deplete first in the region where photon decelera- 
tion is the largest, whereas a laser with the optimal shape 
loses all its energy in a depletion length due to uniform 
photon deceleration, thus enhancing instantaneous en- 
ergy deposition and wakefield. For a numerical example, 
we consider the optimal and Gaussian pulses of total en- 
ergy 0.5J, wavelength 1/im and cylindrical radius 10/im 
in a plasma with Up — lO^^cm^'^. The transformer ra- 
tio, the maximum wakefield, the required pulse length. 



and the corresponding peak oq are shown in Fig. Q as 
functions of depletion length. 

From Fig. || we see that the transformer ratio and the 
maximum wakefield are consistently larger for shaped 
pulses. In fact, the lines for optimal pulse wakefield 
and transformer ratio represent theoretical upper lim- 
its for all pulses of given energy. The Gaussian pulse 
achieves a maximum transformer ratio when its length 
(measured here as FWHM) equals 1/2 of the relativis- 
tic plasma wavelength. The effects of shaping are espe- 
cially prominent for longer pulses, where Gaussian pulse 
yields almost no wake excitation due to plasma oscilla- 
tions inside the pulse that cause part of the laser photons 
to absorb energy from the wake. On the other hand, a 
shaped laser postpones plasma oscillation until the end of 
the pulse, and all photons decelerate uniformly. For very 
short pulses, the differences between the two shapes are 
minimal. This is due to the fact that very short Gaus- 
sian pulses of fixed energy asymptotically approach the 
delta-function limit of the short optimal shape. 
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FIG. 2. Comparison of the transformer ratio, maximum 
wakefield, pulse length, and maximum normalized vector po- 
tential in shaped (diamonds) and Gaussian (triangles) pulses 
of equal depletion lengths and constant pulse energy of 0.5J. 

Although short pulses generally produce the largest 
wakefields, their efficiency is close to minimal possible, 
as the depletion length decreases faster than increase in 
the wake. Therefore, the choice of the appropriate pulse 
shape for LWFA stage will depend on specific experimen- 
tal conditions. If the laser-plasma interaction distance is 
limited by instabilities, diffraction or dephasing, then in 
order to maximize the electron energy gain one should 
try to achieve the largest accelerating gradient, which 
can be accomplished with ultrashort pulses. If, how- 
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ever, the interaction length is less constrained, such as 
the case for propagation in plasma channels |13|, then 
using a finite-length shaped pulse will result in a greatly 
improved overall energy gain per stage. An added ben- 
efit of pulse shaping is the suppression of modulational 
instability that affects unshaped pulses that are longer 
than plasma wavelength. When all photons redshift, or 
"slow down", at the same rate, different laser slices do 
not overrun each other, and the ID laser self-modulation 
is suppressed. 




FIG. 3. Laser intensity (shaded) and the associated photon 
deceleration ( — for pulses of the same total energy and 
characteristic depletion length in the order of decreasing a. 

As the optimal pulse shape is associated with a delta- 
function precursor, the feasibility of such a structure may 
be a concern. We note that the purpose of this precursor 
is to bring the photon deceleration from zero in the qui- 
escent plasma before the laser to a finite value at the 
beginning of the main pulse. This can also be achieved 
with a more physical prepulse, whose shape can be found 
once a smooth function is chosen. For our example 
we choose k{C^ = ko[1 -I- tanha(C — Co)]/2, where a is a 
steepness parameter and is an arbitrary offset. The 
corresponding laser shape is then: 
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where Ci = C - Co < C/ and x(Ci) = 1 + (^o/Za) In ^ - 
Ko[Ci + ^ lii(coshQ;Ci)]/2. As before, the pulse length C/ 
can be found from the available pulse energy. For a step- 
function photon deceleration (a — >• oo) expression ( |ll| ) 
asymptotes to equation (|^). However, for finite values 
of a the delta-function precursor spreads out and can 
even disappear as shown in Fig. ^. The family of shapes 
given by (11) is better suited for the finite-bandwidth 
laser systems that have a lower limit on achievable fea- 
ture size. The values of the maximum wakefield for pulses 
in Fig. ^ are within few percent of the value for an equiv- 
alent delta-function optimal pulse because the bulk of 
the modified laser pulse still experiences constant maxi- 
mal photon deceleration. The wakefield further degrades 



with longer rise times of k(C)- It is also possible to con- 
struct optimal shapes that propagate in a pre-existing 
plasma oscillation and act as wakefield amplifiers [ pT| . 
Such shapes also do not require delta-function precur- 
sors. 

Several issues should be addressed before the laser 
pulse shaping concept can be fully utilized. Even without 
the delta- function precursor, the finite laser bandwidth 
will necessarily smooth out sharp rises and falls of the 
optimal pulse shape. Although we do not anticipate ad- 
verse effects when the feature size is much smaller than 
the plasma wavelength, the ID self-consistent laser evo- 
lution and stability of realistic optimal shapes are cur- 
rently under investigation. Another consideration is the 
infiucnce of the laser-plasma interaction in the trans- 
verse dimension on the evolution of the pulse. Many of 
the laser-plasma instabilities are seeded by the wakefield- 
induced perturbations of the index of refraction. As we 
have demonstrated in this Letter, the nonlinear index 
of refraction can be effectively controlled through laser 
shaping, thus suggesting the method of delaying the on- 
set of these instabilities. Whether this approach increases 
the growth rates of other instabilities, particularly in the 
transverse dimension, remains to be investigated. 
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